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Abstract 

Statistical properties of the temporal distribution of polarity reversals of the geomagnetic field are 
commonly assumed to be a realization of a renewal Poisson process with a variable rate. However, 
it has been recently shown that the polarity reversals strongly depart from a local Poisson statistics, 
because of temporal clustering. Such clustering arises from the presence of long-range correlations 
in the underlying dynamo process. Recently achieved laboratory dynamo also shows reversals. It 
is shown here that laboratory and paleomagnetic data are both characterized by the presence of 
long-range correlations. 

PACS numbers: 
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I. INTRODUCTION 



The observation of the paleomagnetic data l|-[3| have shown that, unhke the solar mag- 
netic field, where the polarity reversals are strictly periodic, geomagnetic measurements 
of the last 160 million years present rather sudden and occasional polarity reversals. The 
reversal process is normally very rapid with respect to the typical time interval between 
successive reversals, which may range from 10^ up to 10^ years . Recent works on 

data analysis, experimental dynamo and theoretical modeling have inproved the knowledge 
of the Earth dynamo. However, the main fundamental questions concerning the polarity 



reversals still remain unanswered ^^\. The nature of the triggers (external or internal 
to Earth) and the physical mechanisms giving rise to the reversals, the reason for the long 
time variations in the average reversal rate (cf. e.g. {2, 8|), are still open problems. 

The sequence of geomagnetic reversals (see the example from the CK95 database jsl shown 
in Fig. [1]) seems to result from a of a stochastic process. The same behaviour is observed for 
experimental dynamo and from numerical simulations [7]. While experimental dynamo 
is a recent excellent achievement, the numerical approach, namely the direct solution of the 
Maghetohydrodynamics (MHD) equations (see 



1 appn 



is still far from being satisfactory 



for a statistical analysis. However, reversal 



fied models, such as few modes models 
two metastable states 



12- 



15 



s are also observed in field resulting from simpli- 
14| . models of noise- induced switchings between 



17| , or mean- field dynamo models with a noise-perturbed a profile 



Recently, it has been shown through a simple statistical analysis, that the reversals of 



the paleomagnetic field are not random 



18N20l|. namely the statistics of interevent times 



(At = tj+i — ti, where ti is the time of the i-th event of the record) departs from a Poissonian 
distribution (namely an exponential law -P(At) = Xexp{—XA(t)), where A represents the 



reversal occurrence rate 



fly 



211), including a non- stationary Poisson process, in which 



case a power-law distribution could arise from the superposition of Poisson distributions with 
time variable rates A(t), see This result shows that geomagnetic reversals are clustered 
in time, probably because of presence of memory in the process generating polarity reversals. 

Here we show that experimental dynamo reversals also are characterized by correlations 
and clustering, suggesting that the reversal process is a universal property of dynamo, which 
does not need any external triggering. 



3 




FIG. 1: Bottom: Polarity of the earth's magnetic field (from today) as in the CK95 record (partial). 
The black bars are the normal (present) polarity. Top: the probability density function P{At) of 
persistence times At for CK95 database (statistical errors are shown as vertical bars). 

II. LOCAL POISSON HYPOTHESIS AND PALEOMAGNETIC DATA 



In this section we will describe the statistical tool used in this work to test, as a zero- 
th order hypothesis Hq, whether the observed sequence is consistent with a Local Poisson 



Process. The reversals rate profile A(t) being in princip 
independent on it. A method introduce d in cosmology 



e unknown, the test should be 



solar flares 



24 



23| and more recently used for 



25| geomagnetic activity 26)| . random lasers in liquid crystals 27|], and 



stock market analysis 28| will be used here. Consider the time sequence of reversals as a 
point-like process, and suppose that each reversal occurs at a discrete time tj. The suitably 
normalized local time interval h between reversals can be defined by introducing 5ti as 



5ti = min{ti+i - tf, U - ti_i} 



(1) 
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and Tj by 



ti-i — ti-2 if Sti — ti — ti-i 



6ti and Tj are then the two persistence times following or preceeding a given reversal at ti. 
If the local Poisson hypothesis Hq holds, both Sti and Ti arc independently distributed accord- 
ing to an exponential probability density: p{dt) — 2Ai exp(— 2Aj5t) and p(r) = Ajexp(— Ajt) 
with local rate Aj. The distribution of the variable h defined by 

= 2^1^ 

will not depend on Aj. 

For the surviving function of the probability density 

poo poo ng{x,H) 

P{h >H)= P{h)dh = / dx2Xe-^^'' / dyXe-^y (4) 
Jh Jo Jo 

where P{h) is the probability density function of h and 



1-1 

H 



g{x, H) — 2x 

it can be easily shown that, under the hypothesis Hq, 

P{h>H)^l-H, 

that is, h is a stochastic variable uniformly distributed in /i e [0; 1]. 

In a process where r^'s are systematically smaller than 25tiS, clusters are present and the 
average value of h is greater than 1/2. On the contrary, when the process is characterized 
by voids, the average value of h is less than 1/2. From time series, it is easy to calculate the 
surviving function P{h > H) and the probability density function P{h). 

The test described above has been recently applied to four different datasets of geo- 
magnetic polarity reversals, including the already mentioned CK95 citeprl,pepi,gafd. The 
probability density function P{h) is reported in Fig. reffig2 for the CK95 datasets. A sig- 
nificant deviation from the uniform distribution was observed in all the datasets, due the 
presence of clusters. 
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FIG. 2: Probability densities P{h) of the stochastic variable h and corresponding surviving func- 
tions P{h > H) for all the empirical datasets. The theoretical probability expected under a Poisson 
statistics is also shown. 



III. EXPERIMENTAL DYNAMO 

The dynamo laboratory model Qj mimics an alpha-omega cycle where part of the dynamo 
cycle is generated by an external feed-back but the flow turbulence is still included and has a 
leading role. In order to achieve this in a simple laboratory dynamo, we relax the requirement 
that the current path be fully homogeneous, and we effectively prescribe an alpha mechanism 
by which a toroidal magnetic field generates an induced poloidal one. However, the omega 
poloidal to toroidal conversion still results from a fully turbulent process. Our experimental 
fluid turbulent dynamo is very much inspired by a variation of the solid rotor dynamo 
proposed by Sir Edward Bullard in the early 20th century, and described in figure 1. Two 
coaxial disks counter rotate at a rate Q. When an axial magnetic field Bq^ is externally 
applied, the flow differential rotation induces a toroidal fleld Bg, this is the omega effect. 
The value of this fleld is then used to drive a linear current amplifler in the loop that 
generates Bq^. The poloidal to toroidal conversion is entirely due to the fluid motion, and 
incorporates all turbulence effects. It has been extensively studied in previous "open loop" 
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FIG. 3: (a) Omega effect : the differential rotation on a von Karman flow advects and stretches an 
externally applied axial field Bqz so as to generate a toroidal component Bq. (b) Postive feed-back 
: the amplitude of Bq is used to drive a power source which generates the current in the external 
loop. Two Helmoltz coils are set on either end of the cylindrical flow vessel; Bg is measured in the 
mid-plane by a Hall probe connected to a Bell gaussmeter. The measured value is fed into a linear 
amplifier whose output drives a Kepco current source. In order to explore the role of the turbulent 
fluctuations, the amplifier has separate channels for the DC and flcutuating parts of the induction. 



experiments. When Bqz is externally fixed, one has Bg = kRmBoz where Rm = R^^l/X is 
the magnetic Reynolds number (with A the magnetic diffusivity of liquid Gallium) and k is 
a "geometric" constant which in our experiment has been measured of the order of 0.1. The 
toroidal to poloidal conversion is then obtained by feeding the axial coils with an electrical 
current linearly driven by a signal proportional Big, so that B^^ = aGBg which reinforces 
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Boz with G an adjustable gain. In such a closed loop setup, one then has — aGkRmBoz, 
and self-sustained dynamo is reached as Q > 0,'^ — X/GkR"^. Clearly, the adjustable gain of 
the linear amplifier allows to adjust the value of fl'^ to an experimentally accessible range. 
At this point it should be emphasized that although the feed-back scheme is very similar 
for the BuUard rotor dynamo and for our fluid experiment, the expected dynamics is much 
richer because of the strong fluctuations in the turbulent flow, where Reynolds numbers in 
excess of 10^ are reached. Indeed, the von Karman flow is known for its complex dynamics, 
presenting not only small scale turbulent fluctuations but also large scale ones - for instance 
fluctuation up to 114% for the differenteial rotation effect has been reported). Compared 
to the 1963 pioneering experiment of Lowes and Wilkison with solid rotor motions, the 
study here fully incorporates fluid turbulence and the associated fluctuations of magnetic 
induction. The role of these fluctuations, inherent to large Reynolds number flows, remains 
one of the mysterties of natural dynamos, and of noisy instabilities in a braoder framework. 

In this experiment the value of the magnetic fleld at saturation Bs^t is flxed by the max- 
imum current that can be drawn from the power amplifler driving the coils. We measure 
Bsat ~ 30 G, a value such that the Lorentz forces cannot modify the hydrodynamic flow 
-since it yields an interaction parameter of the order of 10^^. The saturation of the insta- 
bility is therefore driven by the amplifler non-linearities rather than by the back-reaction of 
Lorentz forces on the dynamical velocity fleld. As a consequence, the B^ component of the 
generated magnetic fleld saturates at the same mean amplitude Bgat for all rotation rates 
{Bsat corresponds to the magnetic fleld generated by the coils when the current source is 
saturated), the saturation amplitude of the toroidal fleld S^sat = kRmBsat linearly increases 
with Vl. 

Another noteworthy observation is that the presence of turbulent fluctuations plays a 
crucial role in the triggering of the magnetic fleld reversals. In the experimental results 
reported here, the current source is driven by an amplifler whose input is Bq -\- g.b'g, with Bg 
the low pass DC component of Bq and b'^ its AC fluctuating part. This arrangement allows 
to study separately the role of slow variations and turbulent fluctuations in the feed-back 
loop. In the results reported in this article, we have set g = 1.18. A homopolar dynamo, i.e. 
without reversals, was obtained for smaller values of g or when the b'g input in the amplifler 
was replaced by a synthetic gaussian white noise (even with a high amplitude) . 

We show here the results of the /i-test obtained in a reahzation (serie27) with Q — 12 Hz, 
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FIG. 4: Probability densities P{h) of the stochastic variable h and corresponding surviving func- 
tions P{h > H) for the experimental dataset described in the text. The theoretical probability 
expected under a Poisson statistics is also shown. 

and cutoff frequency fc = 600 niHz. Similar resutls were observed with different parameters, 
and this study is left for more extended work. Figure reffigS shows the reversals surviving 
function in the case described here. The behaviour is very similar to the paleomagnetic case, 
indicating again presence of clustering and correlations, rather than a random behaviour. 
This indicates that the mechanism responsible for the clustering is present in both dynamoes, 
suggesting some sort of universality of the process. 



IV. CONCLUSION 



In this short paper, the statistical properties of persistence times between geomagnetic 
reversals have been investigated. We performed a statistical test which showed that ge- 
omagnetic reversals are produced by an underlying process that is far from being locally 

I I 

Poissonian, as recently conjectured by [22]. Thus, the sequence of geomagnetic reversals is 
characterized by time correlations. As spontaneous reversals of the geodynamo field have 



been observed in high resolution numerical simulations 10 



lo|,liil. 



the main results contained 
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in this paper seem to indicate that such reversals could be related to the non-linear nature 
of the turbulent dynamo. In order to confirm this conjecture, we performed the statistical 
test mentioned above on recent results from laboratory dynamo. Our analysis has shown 
that the departure from Poisson statistics found in the paleomagnetic data, related with 
the long range correlations introduced by the chaotic dynamic of the system citepepi,gafd, 
are also present in the laboratory dynamo. Such correlations can be associated with the 
presence of some degree of memory in the underlying dynamo process 0, S] which gives 
rise to clustering of reversals. 
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